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Abstract 

The silicon pixel vertex detector is a key element of the BTeV spectrometer. Sensors 
bump-bonded to prototype front-end devices were tested in a high energy pion 
beam at Fermilab. The spatial resolution and occupancies as a function of the pion 
incident angle were measured for various sensor-readout combinations. The data 
are compared with predictions from our Monte Carlo simulation and very good 
agreement is found. 
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1 Introduction 



BTeV is an experiment expected to run in the new Tevatron CO interaction 
region at Fermilab in « 2006. It is designed to perform precision studies of b 
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and c quark decays, with particular emphasis on mixing, CP violation, rare 
and forbidden decays [1]. This experiment takes advantage of two important 
features of the "forward" region: the correlation in the direction of the pro- 
duced b and b, that improves the flavor tagging efficiency, and the boost that 
is exploited in our trigger algorithm based upon the identification of detached 
charm and beauty decay vertices. The unique feature of BTeV is that this al- 
gorithm is implemented in the first trigger level [2]. Consequently, the vertex 
detector must have a fast readout, superior pattern recognition power, small 
track extrapolation errors, and good performance even after high radiation 
dose. Silicon pixel sensors were chosen because they provide very accurate 
space point information and have intrinsically low noise and high radiation 
hardness. 

In this paper we report the results of the 1999-2000 BTeV silicon pixel detec- 
tor beam test and we compare them to Monte Carlo predictions. The beam 
test was carried out at Fermilab in a 227 GeV/c pion test beam. The pixel 
detectors tested were hybrid assemblies of several combinations of pixel read- 
out chip prototypes developed at Fermilab, and single-chip sensor prototypes. 
The main goal of our studies was to measure the spatial resolution attain- 
able along the short pixel dimension for different sensor technologies and for 
different readout electronics configurations. In particular, we have performed 
an extensive investigation of the effects of varying the discriminator threshold 
and the front-end device digitization precision. 



2 Experimental setup 

The data were collected at the MTest beam line located in the Meson Area 
at Fermilab. Fig. 1 shows the experimental set-up. The pixel devices were 
located between two stations of silicon microstrip detectors (SSD's) that pro- 
vided tracking information to an accuracy of about 2 jum in the x direction, 
corresponding to the "small pixel dimension" (50 /im). The pixel hybrid de- 
vices were mounted on printed circuits boards, held inside an aluminum box, 
where their location was determined by precision machined slots. One of the 
pixel devices in the telescope could be positioned in slots at various angles 
with respect to the beam direction. This enabled us to measure the properties 
of various pixel prototypes as a function of the pion incident angle. 

The pixel detectors tested are all from the "ATLAS prototype submission" 
[3], and all have 50 /im x 400 /im pixels. They use the n + /n/p + technology, 
with the pixel electrodes located on the ohmic side of the device. In order to 
achieve good inter-pixel insulation, two approaches have been tried. The u p- 
stop" technique uses p + implants between the pixel cells, deposited through a 
mask with the chosen implant geometry, whereas in the "p-spray" technique, 
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Fig. 1. Schematic drawing of the silicon telescope and data acquisition system in- 
stalled in the thermo-controlled hut. 





Fig. 2. The photograph on the left shows CiS "ST1" bonded to FPIXO. The sensor 
is on top in the photograph. Only pixels in the lower right corner of the sensor are 
bonded to the FPIXO, which is smaller than the sensor. The photograph on the 
right shows SII "ST1" bonded to FPIX1. All sensor channels are bonded to readout 
pixels. 



pixel implants are deposited after a shallow p + layer is deposited uniformly 
throughout the active area. Two of the sensors were produced by CiS, Ger- 
many, and the other three by SII, Japan. Fig. 2 on the left shows an FPIXO 
bonded to the CiS p-stop sensor. The instrumented portion of the sensor is 11 
columns x 64 rows. The CiS sensors (one p-stop ST1 and one p-spray ST2) 
were indium bump bonded to FPIXO readout chips by Boeing North Amer- 
ica, Inc. The SII sensors (two p-stop STl's and one p-spray ST2) were indium 
bump bonded to FPIX1 readout chips by Advanced Interconnect Technology 
Ltd. (AIT, Hong Kong). The depletion voltage, measured through the depen- 
dence of the leakage current on the reverse bias applied, is 85 V for the CiS 
sensors and 45V for the SII sensors. The sensor thickness was about 300/im. 
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The front-end device to be coupled to the pixel sensors must satisfy several 
challenging requirements. In order to obtain the optimal spatial resolution 
achievable with the chosen pitch, analog information is needed. On the other 
hand, it is necessary to transfer the hit information very quickly to the trigger 
processor in order to be able to use it in the first trigger level algorithm. The 
development effort towards the final chip [5] has proceeded along several steps 
of increasing complexity. 

The first iteration, FPIXO, has 12 columns of 64 rows. Each FPIXO readout 
pixel contains an amplifier, a comparator, and a peak sensing circuit. When 
any comparator fires, a FAST-OR signal is asserted. FPIXO provides a zero- 
suppressed readout of hit pixels. The information read out consists of hit row 
and column numbers, together with a voltage level which is proportional to the 
peak pulse height. For these measurements, the analog output was digitized 
by an external 8-bit flash ADC. 

FPIX1 is the second-generation pixel readout chip developed at Fermilab, and 
is the first implementation of a high speed readout architecture. Each FPIX1 
cell includes 4 comparators: one comparator is used to provide sparse readout 
and the other three implement a 2 bit flash ADC for each cell. 

Two scintillation counters located upstream and downstream of the telescope 
determined the trigger, through a coincidence with the FAST-OR signal from 
the FPIXO p-stop hybrid detector. The apparatus was located inside a tem- 
perature controlled enclosure, maintained at rs 20°C. 



3 Pixel charge calibration 

A relative calibration of the ADC response of each cell on a readout chip has 
been performed by injecting charge into individual pixels, sending a voltage 
pulse to a calibration capacitor in each front-end channel. With this method 
the gain and equivalent noise charge are determined up to a scale factor asso- 
ciated with the value of the input capacitor. In order to perform an absolute 
calibration, we have used two x-ray sources (Tb and Ag foils excited by an 
Am a emitter). The X-rays produce known signals in the sensors and lead to 
an absolute determination of the electronic gain as well as the equivalent noise 
charge and discriminator thresholds. 

In Fig. 3, the pedestal subtracted and gain equalized spectrum of a Tb X- 
ray source measured for an FPIXO-instrumented detector is shown. The single 
channel calibration curve fits the data points very well. FPIXO contains two 
different input cells characterized by different gains ( "high gain" and "stan- 
dard gain"). In the following front-end characterization, the quoted central 
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Fig. 3. Calibration for the FPIXO p-spray (CiS) sensor, (a) single channel calibration 
curve with a nonlinear fit superimposed; (b) Tb x-ray spectrum. The arrow indicates 
the K a peak. 

value gives the average quantity and the uncertainty gives the distribution 
rms within a chip. For most of the data taking, the discriminator threshold 
for the FPIXO p-stop was set to a voltage equivalent to 2500±400 e~ for the 
standard gain cells, and 1500±230 e~ for the high gain cells. For the FPIXO 
p-spray device the corresponding thresholds were typically 2200±350 e~ and 
1250±160 e~~. The equivalent noise charge is 105±15 e~ for standard gain 
cells, and 83±15 e~ for for high gain cells of the FPIXO p-stop hybrid pixel 
devices. The corresponding values for the FPIXO p-spray devices are 80±10 
e~ for standard gain cells, and 67±8 e~ for high gain cells. An additional con- 
tribution due to the external buffer amplifier and ADC of about 400±150 e~ 
for FPIXO standard gain cells and 205±95 e~ for FPIXO high gain cells was 
present. 

The absolute calibration of the FPIX1 readout chip response was determined 
by measuring the differential counting rate due to the same two X-ray sources. 
In Fig. 4 the differential counting rate obtained by sweeping the external 
voltage V t ho for a hybrid detector exposed to a Tb source is shown. The 
peak of the derivative of the efficiency curve defines the threshold voltage 
V t h corresponding to the known K a line. The procedure is repeated with a 
Ag source with a line producing a pulse height within the linearity range 
of FPIX1 to make an absolute measurement of the slope of the flash ADC 
response. The threshold used for FPIX1 p-stop hybrid detector is 3,800 e~. 
In addition, the three thresholds determining the digitization bin sizes of the 
flash ADC were 4,500 e~, 10,300 e~, and 14,700 e~. The threshold dispersion 
for each discriminator was fitted to a Gaussian distribution. The standard 
deviation was ~ 380 e~ in each case. The front-end equivalent noise charge 
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Fig. 4. Calibration for the FPIX1 p-stop (SII) sensor: (a) Tb x-ray differential 
spectrum versus the external voltage Vtho- The arrow indicates the K a peak; (b) 
threshold calibration curve obtained from the Tb and Ag K a peaks and the four 
FPIX1 thresholds in nominal operating conditions. 

is 110±30 e~, where the quoted uncertainty reflects the rms spread within a 
chip. 



4 Monte Carlo simulation of the pixel sensor performance 



In order to identify the sensor and readout properties that are optimal for 
our experiment, we have developed a Monte Carlo simulation [4] including 
detail modeling of all the main physical processes affecting the development 
and collection of the signal in silicon. 

Electrons and holes are produced by the energy deposited in the sensor by 
the charged track, and are simulated including excitations, low energy ioniza- 
tion [8], and energetic knock-on electron (5 ray) emission [9]. The algorithm 
will be discussed referring to u n + /n/p + " pixel sensors, where the relevant 
charge carriers are electrons. The drift motion of the electrons along the E 
direction is described by the current density equation: 

Je = -qpe^eE (1) 



where q is the magnitude of the electron charge, // e is the electron mobility 
and p e is the number of free electrons per unit volume. The speed is related 
to the electric field through the drift mobility /x e . An experimental /i e param- 
eterization [10] is adopted, including non-linear effects at high fields. 

The charge cloud spreads laterally due to diffusion. The parameter character- 
izing the drift in the electric field (fj, e ) is related to the parameter describing 
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the diffusion of the charge cloud (D e ) by the Einstein equation: 
kT 

D e = —fi e (2) 



where D e is the electron diffusion coefficient and kT is the product of the 
Boltzmann constant and the absolute temperature of the silicon. The average 
square deviation with respect to the trajectory of the collected charge without 
diffusion is < Ar 2 >= 2 < D > At. Note that in n + /n/p + sensors, the 
region of highest electric field is the farthest from the signal electrodes. This is 
the region characterized by lower mobility and, consequently, lower diffusion 
coefficient. However the collection time is also influenced by the non-linear 
behavior of the drift velocity, thus the diffusion radius is less sensitive to this 
effect. 

Fig. 5 shows the charge spread for track angles of rad and 0.3 rad with respect 
to the normal to the detector plane. Diffusion determines the shape of the 
charge distribution at small incident angles and allows interpolation between 
pixel centers using charge weighting. At larger angles the charge division is 
linear and is dominated by the cluster broadening produced by the track 
inclination. The spatial resolution is determined by the precision of the charge 
measurement and by the pixel pitch. Our pixel detectors are expected to be 
very low-noise devices, as can be inferred by the data discussed previously. 
In order to take full advantage of this feature, the discriminator threshold 
spread needs to be small, as this spread needs to be added in quadrature to 
the intrinsic noise to determine the overall noise performance of the system. 
As mentioned before, the threshold spread measured in the devices used in 
the test beam was about 380 e~ . Therefore noise and threshold spread figures 
are not limiting factors in the detector performance. 



5 Data analysis and results 

5.1 Analysis Method 

The track position at each plane in the telescope is reconstructed using analog 
charge weighting. Each track is fitted to a straight line using the Kalman filter 
technique and tracks with a good x 2 are retained for further analysis. The 
Kalman fitter includes information from strip detectors and additional pixel 
sensors, but excludes the device under test. The track parameters allow the 
determination of the position (x T , y T ) at which the beam intersects the sensor 
being characterized. Although the pixel detectors measure two coordinates, the 
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Fig. 5. Charge distribution for a track angle of rad (top) and 0.3 rad (bottom) 
with respect to the normal to the detector plane. 

discussion will focus on the resolution that can be achieved in the direction 
with smaller pitch, x. 

In order to predict the track impact point Xt, it is important to align the 
individual detectors in the telescope. We use a right-handed coordinate system 
x, y, z, where z corresponds to the beam direction, y is oriented along the 
vertical direction and x is the horizontal axis. Each plane is defined by 3 
offset parameters 8x, 5y and 5z and three angles a, (3 and <fi defining their 
orientation through three rotations around the x, y and z axes respectively. 
We have developed two alignment procedures, implemented using the first 
1000 events in a run. A "manual" technique finds the optimal values of some 
parameters, while maintaining geometrical parameters that are well measured 
at their known value. In a first iteration, translational offsets are determined 
such that the residual distributions are centered around zero. Next, the two 
relevant pixel plane angles, (5 and (f> are determined by minimizing the residual 
widths in the x and y directions. A second, automatic iterative alignment 
procedure has been developed performing a global minimization of the \ 2 °f 
the fitted tracks. This procedure guarantees that all the runs are aligned in a 
consistent fashion. 

The x coordinate measured in the pixel plane is determined in two steps. First 
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we obtain the "digital" coordinate: 

xr + xl , q x 
x d = ^ , (3) 



where (xr) is the local coordinate of the left-most (right-most) hit in the 
cluster. The "left" pixel is the one with the smallest x coordinate. 

Subsequently we refine this measurement with an empirical correction term 
expressed as a function of the variable rj [6] defined as: 

V = : , (4) 

Qr + Ql 



where q# and q^ are the charges deposited in the pixels located at the right 
and left boundaries of the clusters. Thus the measured position is given by 
x p = xd + f{rj). The correction f{r]) is a function of the number of pixels in 
the cluster and the track angle, and is determined by fitting an independent 
data set. 

The measured spatial resolution of a pixel plane is the difference between 
the reconstructed pixel position x p and the track position xt, known with 
an accuracy of about 2 jum, limited by multiple Coulomb scattering from the 
material in the telescope and the intrinsic hit resolution of the various tracking 
stations. This prediction uncertainty is subtracted in quadrature from the 
measured spatial resolution in the results shown below. 



5.2 Charge signal distributions 



The 8 bit analog information provided by the FPIXO-based readout electronics 
allowed us to measure several features of the charge signal in various sensors. 
The measured pulse height distributions are fitted using a Landau function 
convoluted with a Gaussian [11]: 



(E-E 1 ) 2 



m = N J dE ^T—e ' (5 » 



-oo 



where N is a normalization constant, is the Landau probability distribu- 
tion, with Ao = —0.223 and E mp is the most probable energy loss. For 300 \im 
of silicon we have £ m 5.35 KeV, E mp « 84.02 KeV, and a g w 5.95 KeV. The 
study of the charge collected in single pixel clusters for tracks at normal inci- 
dence determines the straggling function in our sensor. This can be compared 



9 




20 40 60 80 

Pulse Height (x 1000 e) 



100 



Fig. 6. 0° track pulse height distribution for a CiS p-stop sensor bump-bonded to 
an FPIX0 readout chip. The Blunck-Leisegang curve [11] fit is superimposed on the 
plot. 

with model expectations and data from silicon strip sensors. The fit includes 
N/£, E mp + Ao, £, <J g as free parameters. Fig. 6 shows the pulse height distri- 
bution for an FPIX0 p-stop sensor. The full width at half maximum (FWHM) 
is 10,000 e - . 



The agreement with the parameters of the theory and with previous measure- 
ments in strip detectors [12] demonstrates that we are achieving full charge 
collection for p-stop sensors. We can use the mapping of the average pulse 
height as a function of the track position to single out regions with poorer 
charge collection properties. Fig. 7 shows the average pulse height as a func- 
tion of the track position for an FPIX0 p-spray hybrid detector. There is an 
obvious dip at the boundary between two pixels in the y direction, indicating 
a strong charge collection inefficiency A less pronounced dip at the interpixel 
boundary in the x direction is also present. This collection deficit is conjec- 
tured to be induced by a parasitic path to the bias grid, introduced to apply 
reverse bias to all the pixel cells during wafer testing of the devices [13] . An 
improved design of p-spray sensors has been developed, and is expected to 
overcome this limitation [14]. 

Table 1 summarizes the FWHM and the parameters of the Blunck-Leisegang 
function for various incident track angles. 
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Fig. 7. Average pulse height versus track position for the ST2 CiS p-spray sensor 
bump-bonded to an FPIXO readout chip. 
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Table 1 

Parameters of the Blunck-Leisegang function and FWHM of the collected charge 
distribution for FPIXO p-stop detector for various incident track angles. 



5.3 Pixel occupancy 

We performed a thorough study of the "row occupancy" , namely the number 
of x pixels in a cluster as a function of the track incident angle and the sensor 
and electronics operating conditions. Fig. 8a shows the number of rows in 
a cluster as a function of the angle measured for FPIXO sensors in nominal 
operating conditions as well as our Monte Carlo predictions. The agreement 
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Fig. 8. Cluster row multiplicity fractions for various angles of incidence. The plot 
on the left shows data from FPIX0 p-stop (V&j as = —140^ and Qthr = 2500e~), 
the plot on the right from FPIX1 p-stop (y bias = -75V and Q thr = 3700e"). The 
curves are our Monte Carlo predictions. 

between data and simulation is excellent over the whole angular range studied. 
The FPIX1 based hybrid detectors have been studied with the same procedure, 
and the results are shown in Fig. 8b. 

5.4 Position resolution 

We studied the spatial resolution achieved with different sensor-readout elec- 
tronics combinations for six track incident angles (0, 5, 10, 15, 20 and 30 
degrees) . The bias voltage was 140 V for CiS devices and 75 V for SII devices, 
corresponding to overdepleted sensors. The data sample was about 50,000 
events per measurement. We also studied the sensitivity of our results to key 
parameters such as the bias voltage applied to the sensor and the discriminator 
threshold. We collected about 10000 events for each operation condition. 

Only events containing a single reconstructed track, characterized by a clus- 
ter including at least two strips in the most upstream and downstream ad- 
measuring SSD planes are used in this analysis. 

Fig. 9 shows the measured 77 distribution for 10° tracks. An entry is made in 
these histograms only if adjacent pixels in a single column are hit. The curve 
shows Monte Carlo predictions for a track angle offset by 1.8° with respect 
to the value determined with the automatic alignment procedure described 
before. The 10° angle is chosen to determine the global offset reflecting a 
rotation of the beam axis with respect to our apparatus. This is because at this 
angle most of the tracks form 2 pixel clusters and thus the 2 pixel 77 distribution 
has the best statistical power and can be predicted more accurately. The 
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Fig. 9. The curve shows the predicted r/ distributions for 12.6° track angle and 2 
pixel clusters. The measured points are with a FPIXO-CiS p-stop hybrid detector. 

agreement between Monte Carlo and data is excellent, indicating that the 
various factors influencing the charge sharing are well modeled. 

The empirical correction f(j]) is determined by requiring that the average 
value of the residuals of the xd distribution is 0. The average residuals as a 
function of r\ can be fitted with a variety of functions, the simplest one being 
a straight line. In the FPIX1 case, it is convenient to determine a correction 
for each of the 16 possible values of the pair Fig- 10. a shows the 

difference Axp between xt and the digital centroid xo as a function of rj 
for tracks at normal incidence; Fig. 10. b shows the linear fit to the residual 
average distribution f(r]); Fig. 10. c and d show the corresponding plots for at 
10° incident angle. We have determined a set of f(r], n P i X , #t)'s using similar 
plots for each detector. This is done using independent track samples, dividing 
rj into 50 sub-intervals, and for each sub-interval determining the average value 
of Axp. 

The residual distributions for the FPIX0 p-stop detector are shown in Fig. 11. 
Each distribution is fitted to a Gaussian. At small angles the residual distribu- 
tions are non-Gaussian, because single pixel clusters are the dominant multi- 
plicity and have a flat residual distribution. In addition there are non-gaussian 
"tails" due to the emission of 5-rays, that are discussed below. Nonetheless, the 
Gaussian standard deviations provide a commonly used measurement of the 
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Fig. 10. Scatter plots of the digital residual versus the w: (a) for 0° and (c) for 10° 
track angles; (b) and (d) show the corresponding corrections f(n) extracted from a 
linear fit to the average of the digital residual distributions. These data correspond 
to a FPIXO-CiS p-stop hybrid detector. 

spatial resolution, that gives a reasonable parameterization of the dominant 
component of the residual distribution. 



Fig. 12 shows the resolution as a function of angle. Two curves and data points 
are included: the solid line and circles show prediction and measurements done 
with an external 8 bit ADC; the dashed curve and triangular data points 
illustrate the results obtained using only x D , the binary reconstructed position, 
to simulate digital readout. Note the excellent agreement between simulation 
and data. The clear advantage of the analog readout is also evident. The 
binary interpolation is less accurate and features pronounced oscillations in 
the spatial resolution as a function of the track incident angle. 



5.5 Resolution versus bias voltage and readout threshold 



Data were taken with a variety of sensor bias voltages and readout thresholds. 
Fig. 13. a shows the predicted sensitivity of the spatial resolution on bias volt- 
age for a p-stop detector bump bonded to an FPIX0 chip with 8 bit analog 
readout. Fig. 13. b shows the corresponding data. A noticeable improvement 
in the resolution is obtained at small track angle when the reverse bias is close 
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Fig. 11. Residual distributions for the FPIXO p-stop detector. oq is the standard 
deviation of the Gaussian fit to each residual plot. 



to the depletion voltage. This is because the longer collection time with lower 
bias voltage allows more charge sharing, thus reducing the percentage of the 
single pixel clusters. 

The discriminator threshold is one of the front-end electronics parameters that 
has considerable influence on the spatial resolution. Fig. 14 (top) shows the 
predicted effect of increasing the threshold for an incident angle of 300 
mrad, for analog and digital readout. The bottom plot shows the fraction of 
events having N pixels hit for a given threshold. For instance, for a threshold 
of 2000 electrons, about 50% of the events have 3 pixels hit, and about 50% 
have 2 pixels hit. As the digital clustering algorithm exploits the informa- 
tion provided by the number of pixels in a cluster, its accuracy is best when 
there is an almost equal population in two different cluster sizes: one cluster 
size corresponding to a track incident close to the pixel center and the other 
corresponding to incidence close to the boundary between two pixels. In the 
analog readout case, the accuracy of any position reconstruction algorithm is 
degraded as the threshold increases. Note that at increasing thresholds the 
efficiency becomes smaller, as illustrated by the curve N=0, corresponding to 
no pixel firing because all the charge signals are below threshold. Fig. 15 shows 
the corresponding measured data points for FPIXO hybrid detectors. 
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Fig. 12. Position resolution as a function of beam incidence angle for the FPIXO 
instrumented CiS p-stop sensor. The curves represent the predicted resolution Qthres 
= 2500 e™. The oscillating curve is the simulated digital resolution; the lower curve 
assumes 8-bit charge digitization. The solid circles are the Gaussian cr's of the 
residual distribution, and the triangles are the cr's extracted from fits to residual 
distributions made without using charge sharing information. 




Fig. 13. Spatial resolution as a function of bias voltage for (a) Monte Carlo simula- 
tion and (b) measurement. 
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Fig. 14. Influence of threshold on the spatial resolution (top), and cluster size for 
tracks at 0.3 rad incidence angle (bottom). 

5.6 Effective depletion depth 



The effective depletion depth of the detector can be measured from data at 
large track angle. We used two different methods. The first method has been 
originally proposed by ATLAS [16] to measure the effective depletion depth 
of irradiated and non-irradiated sensors. It consists of a track-pixel position 
correlation estimator, where for each pixel over threshold whose center is at 
position Xj the distance Zj between the track and the backplane is calculated 
with the formula: Zi = (a;, — Xi nc ) x tan 9, where x inc is the x coordinate of the 
entrance point extrapolated by the fitted track and 9 is the track incidence 
angle. The z, distribution is flat between approximate and the effective 
depletion depth. The locations of rising and falling edges vary with different 
electronic thresholds. In order to estimate the depletion depth, a detailed MC 
simulation is needed. 

The second method consists of a length-charge correlation estimator using the 
pixel charge within a cluster. For clusters including N > 3 pixels, we can 
estimate an effective depletion depth z D = (1 + q q L+ ^ ) N ™mjf P ■> where q ijR is 
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Fig. 15. Measured spatial resolution for various electonic thresholds and incidence 
angles. 

the charge collected by the pixel at the left (right) edge, q m iddie is the sum of 
the charge collected by the N — 2 middle pixels, and p is the pitch along x. 
The distribution of the estimator z D is peaked around the effective depletion 
depth. 

Fig. 16 shows the data and MC simulation for the CiS p-stop-FPFXO hybrid 
detector. Using j3 = 30°, the effective depth is found to be 293. ± 1.2/im and 
305. ± 1.0/iin for the two methods respectively. The difference between two 
methods (12.3/xm) can be used as an estimate of the uncertainties in the two 
methods. Note that this error corresponds to about 1° uncertainty in the track 
incident angle. So the depletion depth is (300 ± 1 ± 13)/xm. 

5.7 High multiplicity clusters 

A small fraction of events have clusters with a number of hits greater than 
expected, given only the track trajectory and diffusion. Fig. 17 shows the 
measured frequency of high- multiplicity clusters for tracks at normal incidence. 
Our Monte Carlo simulation, which includes energetic 5-rays, reproduces the 
measured cluster multiplicity reasonably well. 
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Fig. 16. The effective depletion depth distribution for the CiS-FPIXO p-stop de- 
tector using a) method 1, and b) method 2 as discussed in the text. The dots are 
measurements and curves are MC simulation. 
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Fig. 17. The cluster multiplicity fraction in term of number of rows for FPIXO-p-stop 
detector at normal incidence. The points are data and shaded histogram is the Monte 
Carlo prediction. 

5. 8 Resolution function shape 

The pixel residual distribution (or resolution function) deviates from a Gaus- 
sian in two ways: the first effect is produced by the non-optimal the charge 
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sharing at almost normal incidence that produces a large fraction of 1 pixel 
clusters; the second is related to large cluster size events due to delta-ray 
emission. These effects have been alluded to before and are now discussed 
quantitatively using a sample of tracks at nominal normal incidence. 

Fig. 18 shows the residual distributions for the FPIXO p-spray detector taken 
with the beam nominally at normal incidence. The plot on the left shows the 
residual for one-pixel clusters. The plot on the right shows the residual for 
clusters of two or more pixels, with a fit with the function F(x): 

F(x) = N 1 F G (x)+N 2 F NG (x), (6) 



where F G is a Gaussian, and F NG is defined as: 



F NG (x) 



r-^ for \x\ < r c 

\r c p l l 

for \x\ > r c 

\x\< I I c 



(7) 



A p i is a normalization constant, r c is the half width of the constant term, 
and 7 is the exponent of the power-law dependence. The non-Gaussian frac- 
tion N 2 /(Ni + N 2 ) accounts for 18% of the total number of entries in the 
distribution. 

The naive expectation of a rectangular shape 50 /im for the residual distribu- 
tion for one pixel clusters is distorted by two effects. The population at the 
"edges" of the rectangle is depleted by diffusion, inducing two pixel clusters 
for track incident near the pixel periphery. On the other hand, a broadening of 
this ideal distribution is induced by threshold dispersion and electronic noise, 
as well as the small fraction of mismeasured extrapolated position xt- 

The second factor that makes the pixel residual distributions non-Gaussian 
is 5-ray emission. Low energy 5-rays which stop in one of the pixels crossed 
by the particle skew the charge sharing and degrade the resolution. Higher 
energy 5-rays cross one or more pixel boundaries and distort the position 
measurement even more. We have studied this effect with the Monte Carlo 
simulation described before. Fig. 19 shows a comparison between predictions 
and data. The distributions are shown using a log scale to show the tails of 
the residual distribution more clearly. The simulation accounts for about 1/2 
of the broad component of the residual distribution. This may be in part due 
to instrumental effects such as the charge losses near the pixel boundaries of 
the p-spray devices. 
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Fig. 18. Residual distributions for FPIXO p-spray detector at zero degree. The plot 
on the left shows the distribution for cluster size 1 and the plot on the right for 
clusters with multiplicity higher than 1. 
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Fig. 19. Residual distributions for FPIXO p-spray detector at zero degree. Points 
represent data for tracks at normal incidence for several different runs and the curve 
represents our simulation predictions. Both distributions include only clusters with 
multiplicity greater than 1. 

5.9 Charge sharing across columns 



Charge- sharing between columns occurs for less than 5% of the tracks in our 
data sample. We collected enough data to study this effect only at 0°, where 
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Fig. 20. The plots show residual distributions in x (short pixel dimension) and y 
(long pixel dimension) for events where the charge-sharing occurs across 2 columns 
but charge is collected in a single row. 

the diffusion is the basic sharing mechanism. 

For p-stop sensors the effective charge sharing region is a rectangular ring 
with approximately uniform thickness along the pixel boundary. Fig. 20 shows 
the residual distributions a x and a y of 1 row and 2 column clusters. The 
measured y p is calculated with the procedure described before, using a linear 
eta correction determined from the data. The spatial resolution along y is 
consistent with the x resolution when charge sharing allows for interpolation 
between the pixel centers. 



6 Conclusions 



An extensive test beam study of several hybrid pixel detectors has demon- 
strated that both the sensors and the front-end electronics chosen for BTeV 
perform according to expectations. The charge collection properties of the sen- 
sors studied are well understood. We have shown that 2-bit analog information 
is satisfactory and, consequently, that the final version of the front-end elec- 
tronics, featuring a 3-bit flash ADC will provide the excellent spatial resolution 
needed to achieve the BTeV physics goals. 
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